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ABSTRACT 


Title of Thesis: Interstellar Material in Globular Clusters 
David Dupuy Heerwagen, Master of Science, 1963 
mi@eols Qirected by: Uco van Wijk 
Assistant Professor 

The chances of survival of interstellar material in globular clusters is 
investigated. The observational praileeei ons of Roberts (1960) and Idlis and 
Nickol'skii (1959) are used. The effects of ultraviolet radiation due to 
Pmememollue stars an the center of the clustervand the effects of the radiataen 
@eelein ceneral are discussed in Chapter III. Ihe reacticn due to coliiervone 
of globular cluster clouds with material in the galactic disk is evaluated 
both from momentum and energy influx considerations, and the probability of 
evectton of clouds from the cluster is estimated in Chapter JV. It is sug- 
Beseed that the most serious uncertainty is introduced by the incompleteness 
factor for hot 9 stars in globular clusters and further observations of these 


not stars are desired. 
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CHAPTER I 
iNTRODUCT ION 
Two recent investigations have stirred again the controversy concerning 
Mempessible presence Of obscuring matter in globular clusters. First, Idlis 
ané Nikol'skii (1959) reported sighting tiny globule-like cbjects in the central 
region of the globular cluster closest to the sun, Messier 4. Roberts (1960) 
investigated obscuring, cloud-like objects which were discovered in 12 out of 
S32 high altitude globular clusters. 

The arguments presented by these investigators that the obscuring material 
is an integral part of the globular clusters seem convincing. They both showed 
that the probability of such vacant or obscured regions appearing as a result 
of statistical fluctuations in stellar distribution of the cluster is extremely 
remote, if not impossible. Roberts also demonstrated that for such obscura- 
(ion to be caused by intervening halo clouds would require a space density of 
uch clouds much higher than in the solar neighborhood - a result again extremely 
unlikely. It was also noted that clouds or obscurations do not appear against 
the background of the elliptical galaxies. 

In this investigation then, we start with the assumption that these objects 
are indeed obscuring medium in the clusters, to be analyzed concerning both 
their reaction with the galactic disk and the possible environment within the 
cluster which affects their formation and behavior. From the investigation 
of Roberts, we sense a possibility that, since not all the clusters show obscur- 
ing matter, there is a selection effect operating to make conditions for forma- 
vion of such objects and their retention and persistence within the clusters 
more favorable in certain cases than in others. This approach then suggests the 


possibility that there may be small but significant differences between individual 


Seer alee 








clusters. However, the globular clusters seem to be a remarkably homogeneous 
sroup of objects with respect to certain properties. Von Hoerner's (1957) study 
of globular clusters showed that their standard deviation in luminosity is less 
than 1 absolute magnitude, that their deviation of diameters is approximately 

2% of the mean diameter, and that their internal structures are very similar. 
This homogeneity suggests that we adopt a standard cluster model as the frame- 
work within which to make our order of magnitude calculations concerning the 
itjeerstellar material, Because of the great amount of observation and study 
devoted to the cluster Messier 5, we have generally adopted this object as our 
model. 

Of the many basic assumptions that must be made in an investigation of 
mis sort, perhaps the most difficult is that concerning the physical state 
and concentration of matter in the galactic disk, lying outside the discrete 
Seeuas and galactic spiral arms. Results from the study of optical, inter-— 
pceilar absorption lines, emission nebulae, absorption and obscuration of 
slarlignt, and from radio astronomy work in both the continuum and with the 
Pecan Line have given a fairly consistent picture of the galactic structure 
in the dense regions. However, a search of the literature reveals that con- 
fusion and controversy presently exist on the state and density of matter 
between the arms. 

In his study of absorption lines in the optical region, Stromgren (1948) 
implies that the non-cloud region is permeated with a HI medium of approxi- 
mately 0.1 H atoms cm7’, Spitzer (1949) suggests at the first symposium on 
cosmical gas dynamics that the interstellar medium is in pressure equilibrium 
except for the relatively small volume of expanding H II regions near hot stars. 
his model would require the diffuse intercloud medium to have a kinetic tempera- 
cure of 10,000 SNe Signifying that it is primarily ionized or HII. A summary 


of his model shows: 
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Relative 


Volume o (cm) P 
Medium of Palen 95% 10 K WG 107)? dynes /em* 
Clouds of HI 5% Tost iott oP ttt 
Expanding y : 
Regions of HII 0.5% 10' K To iott 


However, later Spitzer (1954) states that the “atoms in this tenuous intercloud 
medium could be mostly neutral and still at a relatively high temperature". Dufay (1955) 
says that regions of diffuse hydrogen of iO ReoniRe concentration could be HI regions 
of temperature 1000 °K or less, because of the long times (107 years) required to 
reach equilibrium at this density. A gas at this density cools at a very slow 
rate. However, on page l?5he remarks that much of the intercloud medium could. 
emule Finally, the question of whether the medium is predominantly HI or HII 
is discussed at length in the Third Symposium on Cosmical Gas Dynamics held in 
1958. Spitzer (1958) concludes that "in any case, this intercloud medium is 
Probably at a high temperature, so that its pressure is equal to that in the 
geuser, cooler HI regions. Whether this rarified gas between clouds, and not 
near an early type star, is HI or HII is an open question". 

In all probability, the intercloud hydrogen of the galactic disk is not uni- 
Eemmein concentration, temperature, or degree of ionization. However, an upper 


eS, 


limit to its concentration can be set reasonably as 0.1 cm~ or one order of magni- 
muge Lower than the concentration of 1 H atom em~? assigned to the spiral arms by 
radio 2l-cm line and optical, absorption line studies. Concentrations as low as 
0.0L em? can also be considered reasonable. 

In this investigation, we are concerned with the gross dynamics of the encoun- 
ter process between a dense cloud or globule of the globular cluster and a typical 
eroud or the intercloud medium of the galactic disk. To simplify our calculations 


and to enable us to obtain both order of magnitude answers and limiting parameters, 


we have assumed the intercloud medium to be neutral. In fact, by using the 
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macroscopic approach of transfer of momentum and kinetic energy between large 
regions of gas and ignoring the individual particle or microscopic behavior, 

we reason that the question of a neutral or an ionized condition for the inter- 
cloud medium, interacting with the neutral clouds of the cluster, is not a 
sensitive input to the final estimates. Bates and Griffing (119555) have calcu- 
lated collisional ionization cross sections for the neutral-neutral and neutral- 
ionized interaction of hydrogen. They find, that at the lower energies of approx- 
imately the range of our problem, the ionization cross section of the two reac- 
tions are approximately the same. 

Little of theory or experiment, applying to astrophysical problems of 
these densities, has been done for the case of plasma interactions. Alfvén 
(1960) has considered the case of neutral-ionized collision. He analyzes the 
ease where a non-ionized gas collides with a plasma where the kinetic energy 
per particle of relative motion between the particles of the two media exceeds 
PewtOulgzgatlon energy per particle of the neutral gas. He concludes that the 
neutral gas will be ionized and trapped in the magnetic field and the density 
of the plasma will increase. However, implicit in his treatment, and also in 
the two papers by Kahn (1957, 1958) concerning collision of two ionized gases, 
is the assumption that the two gases have comparable densities. We will find 
that the case which ultimately interests us in this investigation is the colli- 
slon of a dense, neutral gas of relatively small volume with a highly tenuous 
gas of much Larger extent. 

Kahn finds when two rarified masses of non-ionized gas collide at high 
enough velocity that both gases will be ionized after a relatively shallow 
penetration, and that the counter streaming electrons will be stopped in a very 
Short distance. This results in a stationary, oscillating electron gas through 
which protons are streaming in opposite directions. This so-called collective 


instability of the electrons and not the individual collisions between particles 
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is responsible generally for stopping the counter-streaming gases. 

However, Kahn does assume equal density of particles in each stream. For 
“iemerse where & wide difference exists in the concentrations of the two counter- 
streams, a momentum exchange treatment is appropriate to determine the result- 
ing motion of the arrested, now unified conglomerate. 

eo we will restrict our attention in this problem to the question of what 
Peppers LO the cloud of a globular cluster with respect to the cluster's 
standard of rest when it impacts with the medium of the galactic disk. We will 


neglect considerations of individual particle reaction. 





CHAPTER IT 


DISCUSSION OF OBSERVATIONAL DATA AND ASSUMPTIONS 


Idlis and Nikol'skii (1959) detected small obscurations in the central 
regions of M4. These globules appear twice as numerous photographed in the 
blue as they do in the red, thus following a rough 1/ absorption law and indi- 
@2oime sche presence of dust particles characteristic of the plane of the Galaxy. 
They estimated the following average properties for the globules: radius = 0.02 pc, 
optical thickness (t) = 3, mass of dust = 0.001 solar mass (O), total mass = 0.10, 
concentration of hydrogen atoms = 107 cm”, total number of clouds = 140. They 
report these globules as concentrated toward the center of the cluster with maxi- 
mum concentration at about 0.5 pe distance from the center. 

On the other hand, Roberts (1960) has reported the detection of larger, 
more massive objects in 12 out of 32 high latitude globular clusters examined. 
A typical cloud is that reported in M 143 as having a projected image of 1 by 3 pce 
and a total photographic absorption of 5 magnitudes. Roberts estimates the mass 
of dust in this particular cloud as 2.60. Assuming the estimated 1% ratio of 
dust to gas in the Galaxy plane to hold for the globular clusters, he finds the 
total mass to be 260 ©. Actually, for simplicity of calculations, this is as- 
sumed to be the weight of the atomic hydrogen atoms. This would amount to a 
hydrogen atom concentration of the order of Oe ewes enero meee ee typical 
cluster which shows these clouds would have about LOOO © tied up in such objects. 
Meret es does not discuss the distribution or location of his clouds, other than 
Lo show by photograph that they generally appear at the cluster edge where the 
apparent density or luminosity appear to drop sharply. However, from his discus- 
sien of the dimensions of a cloud shown in a photograph of M 13, one can deduce 
that the cloud lies a minimum distance of about 8 to 10 pe from the cluster center. 


The procedure used to determine the 5 magnitude absorption is not discussed by 


Roberts and is probably a very rough estimate. 
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It is clear that the foregoing two investigators are concerned each with 
an entirely different object, differing by about two orders of magnitude in 
Ppeimaecensity and radius. Accordingly, any discussion of the possible origins, 
lifetimes, or evolutions of interstellar material in globular clusters should 
concern itself and attempt to account for both classes of objects. Throughout 
this paper, those objects reported by Idlis and Nikol'skii will be called 
globules while those of Roberts will be termed clouds. 

In making our analysis, we will adopt a detailed model of a typical globu- 
lar cluster. The cluster, Messier 3, has been observed and studied in great 
detail by Sandage (1954, 1957) and others. Accordingly, we use this observa- 
ciorel data as a basis for cur study. Also, we adopt Oort and van Herk's (1959) 
dynamical model (d) of M3 since it corresponds so closely to its observed structure. 

Interstellar material could come from any or all oft three sources: (1) resi- 
dual material left over after original condensation of the cluster, (2) material 
shat has been ejected from evolving cluster stars, (3) material that has been 
wacked up from the halo and galactic disk by the accretion of clouds within 
the cluster. As we see later, sources (1) and (3) contribute a negligible 
amount compared to (2). We adopt the results of Sandage for the evolution of 
M4 which show that approximately 100,000 © have been converted into interstellar 
material over about 101° years lifetime. We also assume the material to have 
the cosmic proportion of elements, as adopted for interstellar material in the 
galactic disk. 

Orbits of Globular Clusters= Most of the rough orbit calculations of globular 
clusters are based on the radial velocity measurements made by Mayall (1946) or 

by Kinman (1959). In this investigation we will use the more recent data furnished 
vy Kinman. He finds that the clusters follow highly eccentric orbits, the average 
eccentricity being 0.8. It is interesting to note that a graph in a recent paper 
by Eggen, Lynden-Bell, and Sandage (1962) also shows an average eccentricity of 


about 0.8 for high velocity stars of the subdwarf variety, which from ultraviolet 
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excess measurements, appear to correspond to the main sequence of globular 
elvugeers. 

Matsunami and collaborators (1959) have computed galactocentric distances 
for 94 of the 118 knowmclusters, using Lohmann's heliocentric distances. We find 
the average galactocentric distance to be 8.5 kiloparsecs (kpc) for their 94 
elusters. Since we know that for highly eccentric orbits, an object is most 
apt to be found near apogalacticon, we make a rough estimate that the average 
apogalacticon distance (Q) is 10 kpc. Using the formula Q = a(l+e), we find 
the semi-major axis to be about 5.5 kpc. Using the formula for the latus rectum 
p= a(1-e), we tind p = ekpc. 

Having found the parameters for the average orbit, we use them to obtain 
the average velocity (relative to the galactic frame of refernece) of a cluster 
as it traverses the galactic disk. For this calculation we assume that the 
average orbital axis makes an angle of 45 degrees with the galactic disk, and 


we use the polar equation for an ellipse, r = A cece (where & 


l+e cos & anomaly ) , 


soiving for r at the two points of intersection. Upon finding r = 4.5 kpc and 
Pomme. we Select r = 5 kpc as a typical intersection point on the galactic 
disk. 

We use the Schmidt model galaxy (1956) for the distribution of mass in the 
Galaxy. We find the following values from his graph for potential (in km/sec“) 
in the (6, 2) plane: at Q = 10 kpc, %@ = 30,000 km/sec; at r = 3 kpc, 
dy = 80,000 km/sec”. Adopting the value V@ = 60 km/sec at 10 kpc and using 
the formula a - & = iy O-, we find V = 550 km/sec. 

We have made the above approximations for a "mean" orbit for a globular 
cluster, in order to arrive at some typical value for the velocity of a cluster 
(relative to the galactic center) as it traverses the disk. In reality the 


orbits are not closed as we will see in Chapter V. However, for the high e orbits 


with which we are concerned, the foregoing simplifications will give an adequate 
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estimate. We intend to uSe this typical value in our analysis of the probable 
reaction between clouds of the cluster and the interstellar medium of the disk 
during transit. We must recognize that any such analysis is by necessity, 
order of magnitude accuracy. We must realize that for each individual transit, 
Pamee the disk clouds rotate with the circular velocity of the disk, the rela- 
tive velocity for collision can vary from about zero to approximately 600 km/sec. 
For purposes of our analysis, we will select 300 km/sec as the mean relative 
Terocitcy. We use 2c xX 10° years as the period for the mean orbit. 

Stability of Clouds anc. Globules= We have adopted Roberts's cloud as being a 
spherical object of 260 solar masses and 2 pe diameter. It is interesting to 
consider the gravitational stability of such an object. For this purpose we 
use Jean's gravitational instability criterion (1928), and find a rough limit 
that a self gravitating object must satisfy in order to be gravitationally 
stable and not break up under the result of its own thermal motions. In apply- 
ing any such criterion, we assume the cloud to be isolated in Space. If there 
were a@ surrounding medium at a finite pressure, then the minimum temperature 
for stability would be larger. The most convenient form of tnis criterion, 
Mwitemecan also be derived from the virial theorem, is: NoNer = 2.5x10°. 1 /ie 
This signifies that a body of mass M (in solar masses) at temperature T must 
have a density of more than Nem? of hydrogen atoms to be stable. For the 

260 © cloud, we find the density of hydrogen atoms to be about 2700 a For 
this cioud to be gravitationally stable, its temperature would have to be below 
LO? te 50 °K, Such a low temperature is certainly feasible in an HI region 
(neutral hydrogen) particularly when appreciable obscuring dust is present, 


and there is a strong possibility that any such object is stable. On the other 


hand, the 0.1 © globules fall short of the critical density. However, Idlis 





has shown that these objects should have a lifetime of approximately ie years 
when considering their thermal diffusion, so the globules are probably long- 


lived objects in their present state. 
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ENVIRONMENT WITHIN THE GLOBULAR CLUSTER 


If, in fact, we are convinced that the clouds and globules do exist in 
certain globular clusters, then the next matter of interest is an analysis of 
the environment within these clusters which encourages the formation and 
existence of such objects. Of principle interest is an estimate of the amount 
radiation from the stars beyond the Lyman limit (912 A), which would ionize the 
hydrogen. Again, we use the observational data available for M3. 

Star counts by Sandage for M3 (1954) have shown that the group of stars brighter 
than absolute visual magnitude +3.5 (1.25 ©) is slightly concentrated toward the 
center of the cluster. It is in this group of more massive stars above the main 
sequence that we expect to find those actively evolving toward the white dwarf 
Stage by ejecting material. We might, therefore, expect to find this matter ini- 
tially with the same distribution as these bright stars. Generally, this would 
mean a Spherical distribution, somewhat concentrated toward the center. More 
specifically, we can give Roberts's estimate that a maximum of about 1000 solar 
masses of interstellar material is available between passages through the galactic 
disk. He makes this estimate for a recent epoch by first finding the number of 
cluster stars (AN) lying above the main sequence for M3, using the luminosity 
function of Sandage (1957). He also computes the time (At) necessary for the star 
to evolve through this region of the color-magnitude diagram. Then after finding 
\e the average mass loss suffered by a star during evolution, (using 1.2 © for 
initial mass and 0.6 © for final white dwarf mass) he estimates the mass loss of 
the stars of the cluster for a half period of its galactic orbit (ST) by the 
meoduct: AN . (OE) : ST, Since more than half of these bright stars lie 
within & pe of the center of our M3 model, we assume as a first approximation that 
we have 700 © of diffuse material spread evenly within this radius. 

This approximation allows us to make use of the concept of "Stromgren Spheres" 


2 Wale 222 


from Stromgren (1939) to show the extent of ionized hydrogen around stars of types 
O05 to AO, assuming that the surrounding material is distributed uniformly and of 


dD 


density 1 H atom cm~, These tabulated radii may then be modified for regions 
of different concentration by multiplying by yre/3, 

To make this estimate for our cluster model, we use the computations of 
Minkowski and Osterbrock (1959) for M3. They found that about 15 hot blue 
stars, of effective spectral type 05 and of considerably fainter absolute magni- 
tude (M = 42,3) than ordinary Population I stars of the same effective temper- 
ature, would be responsible for most of the ionization. The uncertainty of the 
data is illustrated by the fact that only one of these 05 stars has been sighted 
and the rest are inferred from an incompleteness factor computed by Sandage (1954). 
These 15 stars are distributed within 28 parsecs of the center of the cluster. 
Since we are concerned in our calculation with the volume within 8 pe of the 
center, we correct for the lower muimber of stars by noting that for the Oort- 
van Herk model (da), the volume within p parsecs contains about 0.6 of the mass 
of the stars brighter than M = +3.5 contained within 28 pe. So we correct the 
Minkowski-Osterbrock equation for the final total ionized volume, V = ee) 


p Ne 


1.1x10 


ne 


(1948) has given a correction to his tabulated radii of ionized spheres for 


where N number density, oe of H atoms, and use V = eon otromgren 
the dilution effect for stars that lie outside the interstellar medium. Using 
this correction, we find we can neglect the contribution of those blue stars 
outside of 8 pe for our order of magnitude estimate. 

Now, using a total average mass of hydrogen atoms in diffuse matter within 
8 pe as 700 ©, we find the average concentration of H atoms is N = 15 on™, At 
this density, we find that the whole volume within 8 pe would be ionized. This 
result is for a uniformly distributed medium, and does not take into account that 


any material wouid be concentrated toward the center as indicated by the potential 


field in Figure 1. Also there could and probably would be turbulence and density 
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fluctuations. Nonetheless, this particular environment does not appear conducive 
to the formation of HI regions, particularly of the size of 260 ©. Ionized 
hydrogen has a kinetic temperature of approximately 10,000 ae which gives a 
random thermal velocity of about 16 km/sec TigtewOuLd Cause 1 bOrescane tom 
the cluster quickly. 

miorder to form a general impression of the effect of more interstellar 
material within a globular cluster, we will hypothesize that 45,000 © of mater- 
ial is available, which would be the case if the material remains: wndisturbed in 
the cluster for & time interval:eqnal ta trree mean orbital periods, 


About 2,000 © of this would ve within 8 pe. We will again assume it to be 


1.1x10° 3 


(pc7) with N now equal 42 cm~, 
5 


we find that less than one-third of the volume within 8 pc, or 630 pe ~, is now 


Eeveaa Uniformly. Using the relation V = 


ionized. We could visualize, in this particular case, regions of ionized mater- 
ial or emission nebulae surrounding each hot, blue star, separated by larger 
volumes of HI gas. For the case betweer the two extremes just discussed, we 
might picture instead a region of temuous, ionized gas escaping around floating, 
smaller volumes of HI. 

The 2,000 © of interstellar hydrogen would be just ionized, if distri- 
buted uniformly within a volume of 1l pe radius. It is interesting for this 
S2eeeeonecstimate the time interval for ionizing the entire 2,000 ©, assuming 
that all quanta beyond the Lyman Limit (912 A) are fully absorbed by the hydrogen 
avoms. To make this approximation, we must first calculate the radius of the 
sub-Luminous O-5 type star, which is the principal ionizing agent. We use a 
formule. from Russell, Dugan, and Stewart (1938) My= ara ao log R-0. 08, 

Prom Stromgren, we adopt an effective temperature for the 05 star of 79,000 ae 


Solving for R, we get 0.4 solar radii. 
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The energy density at the surface of a star of effective temperature, T, 


Meewetven by the Planck law as; 


. te 
Uv = Cc (Ee WAT-1) 


(4) 


To find the energy density in the frequency range beyond the Lyman limit P 
we integrate the above expression: V4, ' 
OO hv — 2 
uth (% (ety 'VFdv 
No 


a 
hv% 


At these high frequencies, we make use of the Wien approximation: & = ? 


27) S RT 
un om [7 ev? dv 


Integrating, we get: 


a _ hi, r a pT 
u= SE (kTw © oe | +3 (42) + 6(ER)* OE) 


ieweeere 11nd the intensity at the star's surface, we use the relation: 


L= 7s (Yor, a 
= 2% (Rive) (1 +3/a+ ©/pz + Sa?) 
where A= hMApT 
L (2"8%ec) = ANR*I = 1.7 «10 €O%ec 


This is the luminosity beyond the Lyman limit for one 05 star of My = +2.3. 








We can check these figures by giving an approximate B.C. (bolometric correction) 
for an 05 star. From Allen (1955), BC = -42.5 + 10 log T + 29,000/T, and for 
T = 79,000 °K, BC = 6.9. Using the value, Mo, = 4.62 for the sun and 


D3 


bow 5. 86x10 = 1.9x10° ergs/sec. This figure indicates 


sy ergs/sec. , we find L 


05 
that our calculated luminosity | beyond the Lyman limit is of the correct order 
of magnitude. 

Within 11 parsecs, we have approximately 11 (05) stars, so that the total 
energy output absorbed by 2000 © of hydrogen is 1.9 x 10° ergs/sec. Each atom 


of hydrogen absorbs 13.54 ev in order to ionize. In 2000 ©, there are approx- 
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imately 2.4 x 10° hydrogen atoms. This means we require: 


2,4 x 100° 


x leo" eve ie x Oe ergs/ev Sone 10'9 ergs 
to ionize all the material. At the energy output rate of 1.9 x 10°° ergs/sec, 


the ionization Lifetime is: 


Bra 10 "Gras 


ut 
ST az IO years 
1.4 1078 ergs/sec * 3x10" 2e5¢ 4 


Cay: 


Such a short ionization time (this is, of course, a minimum limit because 
of neglect of radiation losses) would suggest that the interstellar material 
is almost immediately ionized as it appears and quickly -Lost from the cluster. 
feweve:, &ll of our estimates have been based on the use of Stromgren spheres 
TOrea froup of stars with the interstellar material uniformly distributed. In 
reality, we would not expect this to be the case. Space densities derived for 
M3 from Sandage's counts for stars brighter than My= +3.5 show that these 
brighter stars have a distribution proportional to RES King (1962) has 
recently given an empirical law for prcjected densities of the inner regions of 
globular clusters substantiating this type of distribution. Assuming that the 
initial space density of the material follows the same law as the brighter 
stars, we conclude that the initial distribution of matter would be highly con-] 
centrated toward the center. If we also assume that the few stars responsible 
Pememosu Of the ionization follow the same RO law, then we can see that the 
central regions provide the likeliest environment for breeding HI clouds. The 
Vermle@esOr ionization varies directly as the total luminosity of the ionizing 
stars and inversely as Ne, Since the total luminosity varies directly as the 
number of ionizing stars, we finally can write that the percentage of ionized 
volume varies as Gx /y? Pe: other words, the outer regions tend to ionize 
more easily. The Pere coite AP ments do not attempt to provide any mechanism 
to account for the actual formation of HI clouds, but merely offer a general 


analysis as to the most favorable regions for such formation. 
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An investigation by Spitzer (1942) shows that for a spherical self- 
gravitating system the interstellar material retained by a potential field would 
tend to concentrate toward the center because of the field and the effects 
of inelastic collisions. This effect would accentuate the initial concentra- 
tion. The fact that Roberts's clouds are found some distance from the center 


seems to pose an anomaly. This point will be elaborated further in Chapter IV. 


It is evident from the foregoing discussion that one of the biggest 
detriments to the formation of HI regions is the presente of hot, blue O stars 
With their high proportion of radiation beyond the Lyman limit. This would 
bugees. that a correlation exists oetween the presence and distribution of 
Oeewers and the formation of HT regions. Any such correlation would not be 
Straightforward, however, because onte dense HI clouds would form, then the 
subsequent appearanze of these sub luminous O stars in the process of evolu- 
tion would not necessarily ionize dense regions. Nonetheless, it would seem 
logical that some correlation between the age of stage of evolution of a 
Seeeter and 42 favorable ervironment for the prodicticn of HI regions might 
Stet This point will be discussed further in Chapter \V. 

We next look briefly at the kinetic equilibrium temperatures we might 
expect to find in the HI clouds of the cluster. A series of papers by Spitzer 
and Savedoff (1948, 1949, 1950, 1955) have concluded that the equilibrium 
temperature of the HI clouds in the galactic plane is nearly independent of 
ats density in the range from about 1 to 10° H atoms eee This conclusion 
acsumes that the gas has sufficient time to reach equilibrium with its asso- 
ciated radiation field. 

Berore analyzing the HI regions in a globular cluster, we first might 
uiscuss some of the more important energy gain and loss mechanisms that are 
involved. These were discussed by Eddington (1930) and later, in more detail 


by Spitzer and Savedoff. Throughout the discussion, we should remember that 
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we are dealing with a medium not in thermodynamic equilibrium, and in a dilute 
radiation field. However, the frequent elastic encounters between the particles 
do tend to establish an equipartition of the kinetic energy of translation, ard 
the resulting approach to a Maxwellian velocity distribution can be defined 

by an effective kinetic temperature. 

We can assume that there is virtually no radiation in HI regions beyond 
the Lyman Limit, 912 A. ‘The primary heating mechanism is the superelastic 
collisions of electrons with those atoms having ionization potentials below 
15.54 volts and which are relatively abundant, i.e. C, Mg, S, Fe, Al. 

Heating by cosmic rays, at least in the solar neighborhood, is probably negli- 
eeete fOr HI regicns of 1 atom on," or denser. Tne energy Loss mechanisms 

are (1) inelastic collisions of electrons with ions (particularly CII and SilII), 
(2) inelastic collisions between H atoms and grains (3) inelastic collisions 
between electrons and He molecules, and H atoms and He molecules. 

While Spitzer and Savedofrf do nct consider time-dependent sources, Kahn 
(1955) estimates that collisions between separate clouds could be even .. 
more effective than photoionization of the atoms as a heating mechanism. By 
their analysis, Spitzer and Savedoff arrive at an equilibrium temperature 
between 30 and 100 “K for HI clouds of density 10 aie or more. At low 
densities between clouds where radiative cooling processes are ineffective, 
temperatures upward of 500 °K are postulated. Kahn shows that cloud temper- 
avures can be expected to drop very rapidly down to 500 oe between collisions, 
and more slowly below 500 °K, ‘These analyses are observationally substantiated 
by the 2l-cm measurements which yield a harmonic mean temperature of about 
100 to 125 °K for the HI clouds of the galactic plane. 

in attempting to apply any such analysis to HI clouds in globular clusters, 
we are at once struck with the great uncertainties involved. Certainly, we can 
generally assume that the chemical composition of the medium in the clusters has 
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a iower proportion of the heavy elements than does gas in the galactic plane. 
‘Gectra of the globular clusters are interpreted to indicate metal to hydrogen 
ratios from 1/10 to 1/100 that of the sun. <A lower percentage composition of 
the elements with low ionization potential would reduce the contribution of one 
heating mechanism for HI regions and perhaps lead to lower equilibrium tempera- 
ure S . 

We can also expect the radiation density inside a globular cluster to differ 
from that of the solar neighborhood. We can estimate this value quite readily. 
We use the following notation in calculating the approximate radiation density 
inside Messier 3: 


u = radiation density (ergs/cm’)., 


I! 


r = distance from center of cluster (ecm). 

@e-= velocity of light (3 x one cm/sec). 

T = average luminosity (ergs/sec) of those cluster stars brighter than Mo = +30 5- 

L = total luminosity of cluster (ergs/sec) 

tr) Se uensity Of Stars in cluster Dbiieniter tnean x 4305 (stars/parsee’) 

ata distance r. 

From the Oort-van Herk data for model (a) of M3, we use the value 3.3 for LAY) and 
190,000 solar mass for the total mass ()) or the cluster. Using the value 
5. 86x10? ergs/sec, for the total luminosity of the sun and neglecting the probably 
small bolometric correction for the cluster, we find: 

L = (3.86 x 10° ergs/sec).(190,000).(3.3) = 2.43x10°7 ergs/sec. 
lor the purpose of determining the radiation density of the cluster, we consider 
that most of the Light comes from those stars brighter than M = +3.5, so that the 
density distribution of these stars can be equated to the distribution of light 
ZL 


eeaces. For a single star: U= 





at any point a distance r from the star. 
nrc 


Now, to account for the contribution of radiation of all the other stars to 


this single point, we integrate over all the cluster, or in this case out to 


alow = 





LOO pe, the visible limit. 


Poweet @ Single point in the _. we have; 


ne 100 QV d 
() us (4 (arrv)dr = fe 
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For the total luminosity of the cluster, we have the following expression: 


(2) | = (wg uiTr?V dr 
Dividing (1) by (2), we a u/ 7 Ff Dale 
ly Witc (1 pea dae 


Qur final expression for radiation density is then; 


_u feeder 
22 = “are fer aac 


Figures (2) and (3) from Oort-van Herk depict the space densities (f) of the 





group of stars brighter than My= +3.5 versus cluster radius (r). These graphs 
are based on star counts for Messier 4 by Sandage and von Zeipel. In order to 
find u, we evaluate equation by carrying out the integrations. We split up the 
interval from 0 to 100 into 3 segments: (1) integrating numerically by Simpson's 
Rule from 0 to 4 psc, (2) represent analytically aay x) from 4 to 7 parsecs, 


and (3) asN = from 7 to LOO parsecs. 


, 
Upon evaluating the two integrals, and using L = 2,43x10°7 ergs/sec, we 
fmraewnat u= 2.9 x Ones” ergs/cm. 


This value can be compared with the density of dilute black body radiation 


as calculated for the solar neighborhood: 5.24 x We ergs/em (Dunham, 1939) 


em (.O( xX Noles ergs/cu’ (Eddington, 1930), for which the black body temperature 


of a particle in space is 3.18 %, 


Sons 





. 4 
Meine the formula u = akt', where a = 7.57/7x10 = eree/em des and k is the 


absorption coefficient, we can solve for the effective temperature of the aust 
in the cloud if we can find an expression for k. 


We know from Roberts that the 260 © cloud absorbs approximately 5 magnitudes 


in the visible light (Aw5x107 em). From Wien's displacement law, X 3, 2 


We can solve for ee for the effective temperature. 
; : =) 
We also know that k varies as Pe so k = 2 eB ; 








max 
“5 2 <5 
1 
i eenweretore have k = es x >» Biving us Us & =e r, substituting 
#1 


our value of u = 2.9xl0 ergs/ cn and solving for T, we find To 34 “XK. 
The observational evidence supplied by Roberts (1959) by his 2l-cm measure- 
ments of ML3 implies an equilibrium temperature of LOO aie or below. Hence, we 


' 
. 


are led to conclude that the temperatures of dense HI clouds in globular clusters 


ere sab least as cool as HI clouds in the galactic plane. 
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CHAPTER IV 


PASSAGE OF GLOBULAR CLUSTERS THROUGH THE GALACTIC DISK 


We now turn to an investigation of the interaction between gas clouds and 
globules in the cluster with the medium lying in the galactic plane. Previous 
investigators have simply stated that any gas present in the cluster will be 
ejected by each passage through the disk. 

The galactic disk consists of a thin layer of gas and dust distributed in 
8 non-uniform fashion. There are individual clouds and spiral arms separated 
by regions of Low density. Densities for each feature are: typical individual 
clouds - 10 ators em”? , spiral arms - 1 ectom cm™>, space between clouds ~ 0.1 to 
0.01 atoms oa As previously argued, w2 assume most of the hydrogen to be 
neutral for purposes of the collision process. The ionized gas in the disk is 
concentrated around the individual 90 and B stars. 

soince we are interested in the processes that tend to act on the neutral 
hydrogen clouds and globules in the clusters, we can see that the collision to 
be Studied is the HL-HI impact. 

As long as the mean free path for particle collision is small compared to 
the dimensions of the cloud or globule, we can treat the gross dynamics of the 
encounter by evaluating the momentum transfer to determine if the clouds and 
globules of the cluster are ejected from inertial force considerations. For 


8 


collisions of neutral H atoms, we can use the geometrical radius (~1O7 


216 


compute the collision cross-section, giving O= 3x10 a (Parker, 1958). For 


em) to 


the globules, N,,=10° em”; for whe cloudce N.c=10° on™. Using the relation, 
Mimean free path)=1/NG", we find » glob= 10 n¢ and pietOr ne Nealaelol loxonsle: 
values considerably smaller than the dimensions of the object involved. 

Kahn (1955) has considered the statistics of size and distribution of clouds 
in the galactic plane by using the observational evidence that: (1) a line of 
Sight of 1 kiloparsec intersects an average of 10 clouds, and (2) clouds occupy 


= ca = 
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about 5% Si tie volume of the salactic disk. He concludes that the average 


eee the 


cloud is 3.75 parsecs radius, with a density of 20 H atoms cm 
number density of clouds in space is aie oni Ths resultgaemecoursce, 1: 
characteristic of the solar neighborhood. We readily see that the mean free 
ath for collision is about LOO parsecs. 
From 21 cm. observations, Rougoor and Oort (1960) have described the general 
galactic disk between 5 and 10 kpc from the galactic center as being 220 pe thick 
3 


with a mean density of 1 H atom cm~. The disk inside 3 kpe seems to be much 


thinner, about 120 pz, and the mean density is given as 0.4 H atoms om, They 
comment that the region between 3 kp® and the central disk at 600 pe appears to 
be virtually empty. 

Peaem OUr discussion in Uhnapters Ii and Il, we note the rollowing iniormma-— 
tion concerning the clouds and glooules: (1) In general, we would expect to 
find HI regions toward the center of tk.e cluster in an isolated cluster. 

(2) From the incomplete information now availatle, clouds seem to appear at the 
edge of the clusters (minimum of 8 to 10 pc) and globules toward the center of 
2 cluster (in the one example, M4). 

Let us assume that a globular cluster contains a 260 © cloud of Roberts's 
description at the center of the cluster. Let us then analyze the possible 
effects of a collision with the galactic disk. In this analysis we will use 
fyescimecrent treatments for the two different types of impact involved. If we 
assume that the cluster traverses a mean distance through the disk of 300 parsecs, 
then at 3500 km/sec (where 1 km/sec —an pereee/ ie. years) we can conclude the 
time of transit 1s about 10° years. However, a cluster cloud would transit one 
ot Kahn's typical, dense clouds in the galactic plane (diameter&7.5 pc) in 
about 10 years. If we compare these two characteristic times with the approxi- 
mate orbital period of an object in the potential field of the cluster, 10! years, 


we see the need for different treatments. 
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Momentum Transrer 





For the case of impact with a discrete cloud in the disk, we consider an 
impulse to be acting on the cluster cloud. We will also assume that the colli- 
plemeaic Completely inelastic, the colliding portion of the disk cloud remaining 
fone cne Cluster cloud. 

For the case of impact with the intercloud medium (0.1 to 0.01 H atoms em?) 
of the disk, we compute the momentum transfer per unit time of the colliding 
mecmun, which is the force acting on the cloud or globule of the cluster. 

We will see that one effect of the collision with & disk cloud can be to 
Kmie@ek a cluster cloud out of the cluster center and into a rectilinear orbit. 
metme Other hand, the force exerted by the medium on a cluster cloud effectively 
shifts the whole orbit in the direction away from the direction of motion. 

mie simplest manner in which to visualize this problem is to analyze iv 
Wien respect to the potential well of the cluster and to consider the motion 
as a one dimensional oscillator. We use the following notation: 

V = velocity of globular cluster relative to the Galaxy (300 km/sec), 

eee oistvance cf cloud or globule from center of cluster, 

me=etadius of cloud or globule of ciuster, 

L = diameter of disk cloud (7.5 pc) 

n = number density of H atoms in galactic disk (cm), 


m= mass of hydrogen atom (1. 67x107=" 


7 gm). 
(V)= mass of cloud in cluster (260 ©). 
@(R) = potential as a function of R in the cluster (xm/sec)*. 
A = potential gradient, 


ie number density of H atoms in globular cluster cloud: 


I = Momentum per unit time or force applied to cloud in cluster resulting 
from transit of intercloud medium of disk. 
E = Total energy per unit mass of cloud in globular cluster. (km/sec)*, (with 


respect to cluster frame of reference). 


ee 














Figure (1) depicts our detailed example of the globular cluster with a 260 O 
cloud at the center of the cluster on its first impact with the disk. 

First, we derive a general formula for computing the force per unit mass 
resulting from momentum transfer per unit time from the intercloud medium to 


the cluster cloud on transit of the disk. We have that: F = nr”. V>.nem,( dynes) 





be 
(1) = - Nem,,( grams) 
(1) F/M) = “2 ov 8 ( om/sec“ ) 


Hewmour particular case, where r = 1 pe and N = 2.7. 10° em a V = 300 i a and 
n=0.1 em™?, equation (1) gives uss; F/M= 2, 2x1077 onfsec” or 2. s(n)” pea 
We can also give for our example, the expression for the force as; 
(2) F(dynes) = Waeo=o enema) 
. OT 
or F = 25x10) dynes: (for nea O.ia 
We have the following relationships between total energy per unit mass (E), work 
| | F.(R,-R, ) 
done on the cluster cloud per unit mass ane Poel a enersy Peis 


mass (0) and kinetic energy per unit mass (507), where ap = velocity of cloud in 
its orbit with respect to the globular cluster: 
oa 
E 5 


and Ey is the total energy at some time a later. As shown in figures (1) and (25 


ae where = is the total energy in an orbit at time t, 
R values in the direction away from the direction of impact are positive. We 
also have the following identity relating our energies: 

th OS ordare (where Wis the velocity of the cloud in its orbit relative to 

the frame of reference of the cluster) 

We can see that the value of F/(Q = wn determines to what extent an orbit is 
displaced to the rear because of the steady force exerted by the intercloud 
medium during about a years. In figure (2) we have depicted the effect of 
force per unit mass exerted by both a 0.1 om”? medium (F/M = 2,5(km/sec) pe) 


and a 0.01 em™? medium (F/M = 0.25). We depict the extreme case for both these 


oO. a3 





media in order for a cluster cloud (of Roberts's description) to be retained 
meeere Cluster for an indefinite period of time. For the 0.1 em” medium, the 
imitans orbit is from <4 pe to +14 pe; for n = 0.01 om”, itels frene-20upe 
ioe PC. 

mMersnould note, however, that for our particular case where the time we: 
transit of the disk is 10° years. and the orbicaly perlod in theseluster > ~10! 
years, that the chance of a cloud escaping due to the reaction of the medium alone 
Heeexvremely remote. Once the cluster emerges from 1ts transit of the disk, 
then the F/M slope is practically zero and, as long as the cloud has a tcetal 
energy less than 13% (km/sec)®, it will not escape. 

However, the reaction caused by collisions with discrete clouds in the 
galactic disk of the order of n = LO om”? is extremely sipniricant, fem eeu 
discussion. To calculate the energy change per unit mass (AE), we use the 


principle of conservation of momentum. For the general case we have for the 


Wess Of material intercepted by the cluster cloud: 
(VY). = ne 1 Te Mh, 
The mass ((Y)) of the cluster cloud is: 


(V2 - b/3ar?. Nem, 


V is the relative velocity of impact (300 km/sec), 
AWis the velocity change imparted to the cluster cloud relative to the cluster 
frame of reference. 
We then have: (¥2..v = (M+M:) AW 
Ves 
1+ Mm) 


This gives the expression AV -. 


(4) Atre 


2 
is 


+ 4rN 


din 
Now we derive a general expression for the energy transfer per unit mass per 


collision. We haves 


5a 





o = 4305 (where WV. is the orbital velocity of the cloud at time t;) 
E. = o+3( Whar) (E,, = energy after impulse) 
ietoesives us the general expression for energy transfer; 
2 
Aw) 
oe ee 
5) ite EE, = U,AV+ 


Substituting the appropriate values into equation (4) for the 260 © cluster 
cloud colliding with the n = 10 em”? disk cloud,we find AV=6 km/sec. From 
equation (5) then, we have AE = 6 (43). 

Now, with the use of figure (1), let us follow through in some detail the 
reactions of a 260 © cloud, initially at rest in the center of a globular cluster, 
as it passes through a typical area of the galactic disk. We adopt 0.1 H 
atoms om”? as the concentration for the intercloud medium. We assume that the 
average cluster passes through 4300 pe of the disk. The mean free path for 
collision with one of the disk clouds (L = 7.5 pe, n = 10 em”? ) is’ 100mece 

Upon gradual entry of the cluster into the 0.1 om”? intercloud medium, we 
Bemeeu une cloud to be shifted about 0.2 pe to the rear of the cluster cencer, 
or to R = +0.2 pe. The cluster cloud is likely to be near this point on the 
eerem upon its first collision with a disk cloud. Using our previously calleg— 
lated value of AV= 6 km/sec, we find from equation (5), where U,= 0, AE = ef) 
We indicate this change in energy on figure (1) by a dotted vertical line to 
point (R, , E,) or (+0.2, 19). At this point we assume that about one third of 
the transit period (~3x10° years) has elapsed. Through this point we draw a line 
with the proper slope to represent F/Ne2.5 (= ve", as previously calculated 
for the force per unit mass exerted by the 0.1 on? medium on the 260 © cloud. 
This line represents the locus of (R,E) points for the rectilinear orbit which 


the cluster cloud now follows. 
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The increment of energy, AE = 18 (=), received by the cluster in the 


first collision has given the cloud an impetus toward positive R values. We 
now estimate the distance traversed along this orbit in a time interval, 
3x07 years, until the next collision. We integrate numerically along the 
path from R= 0.2 pe te an R value that corresponds to a time interval of 
3x10? years. We use distance steps of 0.2 pce: 

R4+n( 0.2) aR 


R ATR) 


We then find that the cloud will be at (41.6, 22) at the time cf second collision. 


+o 


a 


As depicted in figures (1) and (2), we follow the cloud through two more 
collisions. After the second collision and a time interval of 3x10? years » 
we find that a AE = 53 (km//sec)* for the third collision raises the cloud to 
@ very high energy total (125 (xm/sec)*) which means the cloud probably 


escapes. 


Peeteion, it will find itself im an orpit, ‘with the cistance limite7 oO a0ou 
7 parsecs. From figure (2), we see that a cloud after the third collision is 
supposedly left in an orvit extending out to about 60 tc 70 parsecs. However, 
it is easy to see that, even for a very small frictional force which could 


As 107" oe the cloud 


reasonably be expected from the hale medium (107 
would be lost from such a high energy orbit. 

While we have worked out only this one example in any detail, it is easy 
tO judge the possible reactions in other situations. 

For instance, a subsequent pass by our cluster in the example above ywhich 
has had two collisions on the previous pass through the disk OR eens before. 
might find its cloud in an orbit at an energy level of about E = 70. The 
initial velocity, for this case might vary anywhere from about VU. = +12 km/sec 


o 


through O and back to -l2. Minus values for velocity indicate a direction for 


amo) a 





ite cloud in the same direction as the motion of the cluster orbit. Using 
formulae (4) and (5), we find that AE can have the following values for each 
ev; for +12 km/sec, AE = +903 O km/sec, AE = +18; -12 km/sec, AE = -5h. 
Hence, we see that in these 5 cases: the cloud escapes, has its energy raised 

; a ; 2 ; 

about 18 (km/sec)“, and has its energy lowered 54 (km/sec)~ respectively. 
statistically, we would expect the cloud to be near an extreme point of its 
Orbit and thus find that the change in energy for an impulse is small. 

Two other factors are worth brief mention: (1) we.could expect some 
elliptic orbits to ensue,for the reason that a cluster on a@ subsequent pass 
could collide with the disk at an angle to the direction of motion of its rectilinear 
orbits; and (2) we might expect some damping of the orbit during the period the 
cluster is away from the disk, because of accumulation of additional medium 
from evolution of stars and the effects of inelastic collisions. 

To estimate the chance of a particular cluster cloud getting through the 
disk without escaping from the cluster, we use the results from the previous 


discussion and apply the Poisson frequency function: 


¢ 
k =3 
f{k) = : 


ae 
3 
where a is the expectation and f(k) is the probability that x = k. For our 
problem, we take a equal to 3 (we can expect about 3 collisions per pass). 
femeien evaluate the individual probabilities for k = 0, 1, 2 and 3 collisions 
per transit. We sum up these individual probabilities, using 5. f(k) Or vier or 


4 4 


and find that the probability of a cluster cloud surviving one disk passage is 


about s. Of course this figure is for the average assumed conditions of the 


13k and would vary for a different structure in any region. 
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We should note that the foregoing development for momentum transfer 
applies equally well for the Idlis cloud. That this is true is readily 
understood from equations(1) and (4) which show that the forces depend 
on the product rN which are approximately equal for the two cases. 

It is clear that any original gases in the globular cluster are 

10 

Swept out at some time during passages through the disk during 10 years o 
We previously mentioned one other possible source of material for the cluster, 
that of sweeping wo the gases from the halo. spitzer gives a value of 

Ly — h : a4 : “ p = 
5xlO™" 4H atoms em” for the halo. Shklovsky and Pikelner (1958) give 

5 : 
10 “en ~. If, for example, a Roberts cloud swept a path during its orbit 
of approximately 350 kp: length, and retained all swept matter, the total 
Mess accumulated would be only acout 20 ©, even by using N = LOS. We 


can conclude that a "snow plow” accumulation is a relatively small source 


Ot material. 
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mrey Transfer 


In addition to the effects of momentum transfer, we must also consider the 
heating etfects of a collision between clouds at these high energies. From our 
discussion in Chapter III of heating and cooling processes in cosmic clouds, we 
learned that if an HI region becomes ionized, it will have an equilibrium temper ~ 
ature of about 10,000 “K and the root mean square velocity of the protons will 
pewapout 16 km/sec. The Oort-van Herk dynamical model for MS gives a velocity 
of escape of 16.6 km/sec at the center ranging down to 11.4 km/sec at 8 pe and 
sO on, as can be surmised from figure 2. We conclude that when hydrogen is 
femeeo FO & Kinetic temperature of 10,000 ie it will quickly expand and escape 
trom the cluster. 

Bates and Griffing (1953) have calculated the icnization cross section for 
the reaction H(1S)+H(1S) = H(1S)+i°+2, Extrapolating their data, we estimate 


Ole 
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this cross section as 10°°‘em for a velocity of about 300 km/sec. Comparing 


o a © ° a) Oo o » ’ o u 
the ionization cross section (~10 “‘em’) with the collision cross section 


16 
em”), we see that an individual hydrogen atom has a very low probability 


(~10~ 
(~107 +) of ionization per collision. The incoming energetic hydrogen atoms, 
relative to the frame of reference of the dense (~10? 1A) 10° em”) cloud or 
globule of the cluster, will be quickly damped to energies below 13.6 ev as 
a result of both momentum transfer and radiation losses by excited atoms. As 
we will see later, the ionization will be mostly confined to a narrow interface 
between two colliding clouds. 

We should keep in mind that collisional ionization is not a very efficient 
process. Alfvén (1960) has proposed that collisional ionization between a 
plasma and a neutral gas is caused by electrons being accelerated to sufficient 
energies by irregularities in the charge distribution created in the plasma. 
in the discussion following this paper, Spitzer remarks that when ionization is 
produced by energetic electrons, the cross section for excitation is much greater 


than that for ionization and 4 to 5 times as much energy goes into light than 
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into ionization. Naturally, for any ionization mechanism, a certain proportion 
of the input energy will be dissipated as light and infra-red radiation from 
excited atoms and from grains. 
We will analyze in some detail the case of collision between a Roberts 
cloud (M10? cm?) and a typical cloud in the disk (n~10 an. diameter LL}7.5z7°). 
First, we should note that at relative velocities of about 300 km/sec, the 
Mmieeracting particles exceed the velocity of sound and we can expect a compressea 
eee. 1ront between two cclliding clouds. The velocity of sound is given by 
a =\)¥ RT, Wiel = 2 7 on fo 8.3x10'. HOre arg. Ve a = lea km/sec; ror 
T = 10,000 “K, a = 17 km/sec. 
We now calculate the total energy psr sévord transferred to the more dense 
Roberts cloud as it passes through the larger -~Loud of the disk. We have for 
energy /unit a 
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ergs/sec 


: . o | 
munit time = B( ar“ xV)x( mm) xv" = 6.4210 


(where M, 2 ko om? for disk cloud, r = 1 pe for Roberts cloud, V = 300 km/sec, 
24 de ene - -2  -L | 
em). This means that ——— or 0.23 ergs em sec § 1s the 
wr 
amount of energy input per unit time per unit arca wnich must be dissipated at 


mm, = eo (x10 


the interface of the Roberts cloud. We note that this value is comparable to 
the average stellar radiation per unit area in a clusters 


uc Per 10. Sx 10-~ ft 
= jj 


=- 0.22 ergs em™“sec 
We should investigate the possible cooling mechanisms which are able to 
aissipate this energy and maintain the layer at a reasonable equilibrium temper- 
eeure., 
Spitzer (1949) has given the following formula for computing the rate of 
energy Loss caused by the inelastic collisions between hydrogen atoms and dust 
grains: 


a 
ae ( ergs/cm’sec) = 9. 4x10 7°" Mm g)N(H)E?( T-Tg) 
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We take tne following values tor: 


oa 


n( g) - concentration of dust particles = LO 


Go eee era ty < en ee ) 


(where radius of particles 10 

N(H) ~ value for Roberts cloud 10° em 

(the above two values give a mass ratio of dust to gas of approximately 19%) 

T = gas temperature assumed to be 10,000 EKG 

T . temperature of the grains is assumed very small (10-50 aon 

8 

[ieerting these values in Spitzer's formula, we find Mig = S. tsclom wanes om ~ Sec 
We have an energy input by ccllision of 7.23 ergs om--seo » Which even if dis- 
tributed throughout a layer cf several mean free paths for collision 

1/NO = 1/10°.3x1072° = 1? many orders of magnitude greater 
than the cooling rate. 

We assume that the impact layer be ween the twe clouds is completely ionized 
and at a temperature of 19,009 “x, Accordi ingly, we will have a recombination 
of electrons and protons with the accompanying radiation in the visible and 


infra-red as the electrons cascade down tu the ground state. We use the follow] 


ing formula to obtain an estimate of the cooling rate provided by this mechanism: 


(6) L ac, (eres ees = NN VU $O(nV ),.D 
where N_ = proton concentration in the impact layer (om?). 
Ne = electron concentration (com?) 
ify = velocity of electron at 10,000 e 
On = capture cross section to n° vh state for hydrogen. 


aN) - average energy radiated per capture to state n, and in cascade 
to ground state. 
D = thickness of ionized impact layer (cm). 


From page 91 of Allen (1955), we take the values given for the recombina- 


tion cross section to the oom state for hydrogen and multiply each value by 


(nv) , and then sum all products to give ea) = 5.7x10~- “ergs om7= 
aoe cach@ el ue Onder Or 102" om2). ; 


Zee 





For the velocity of the electron at 10,000 %k, we useW= 6.4x10' cm/sec. We 
consider that NSN. and sev the equaulon. For eae equal to the energy input. 
We then solve for np; 

ew) (FS O (nv) = energy input from collision 


ND (6. 4x10)(5.7x107?°) - 0.23 ergs em” “seo” > 


22 «= 
(7) or Wp = 6.3x10~°cm™ 
Having obtained ession for N- 4 
aving obvained an expression for D, we now find N from considerations 
of pressure equilibrivm in the collision process. We note that the collision 
picture can be deacribed in the frame of reference in which the ionized shock 
layer between the two clouds is stationary so that a disk cloud (n = 10 cm?) 
moves with a velocity of about v = 400 km/sec. We can readily ComoiGe s rim ou. 
Of transfer of momentum per unit area at the front by the expression: 
2 7,2 2 2 


” O4 oe 2 
nem «v= 10 Ge Vor te (om). (3x10')"em™sec 


Cane ] 
u 


= inet = gin en “sac 
pDince this value is the force per unit area or pressure exerted on the 
layer from the direction opposed to the direction of motion, we can equate it 
tO Our expression for the pressure inside the layer, 
p= qe (where N = number of protons per ae k = Boltzmann's constant 


(1..38x10716 


ergs/degree), and T=10,000 Korn We use a as the concentration in 
thas expression since the electrons contribute only about 1/43 ene tle (ene wis 

aac =e ~~ 4 =3 
Pieesoure. We find N ~10 cm ~. 

p 
Using our result from equation (7) that Np = 6.3x10--em™ , we then solve 

for the thickness of the shock layer, finding D = Ae ae We note that D 
is larger than the mean free path for collision of H atoms (A = 3x10 om) in 

: } 
the dense cloud and comparable to the mean free path (X =3x10> tem) in the disk 


cloud. We also note that it is extremely thin compared to the diameter of the 


cluster cloud (2 pe). 
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It follows that we have a situation in which the two clouds and their 


interface are in pressure equilibrium. Using the expression for momentum 








2 z 
2 | 2 Von (500) al 
easier per unit area, nmr , re see that: = z . 
a h Ue eee eee N 3 
10 
er (U. oe km/sec (relative to frame of reference where shock layer 


cluster 


is stationary). 
We then see that the distance of travel by the shock front relative to 


50 


She Cluster cloud during the time of collision is equal to 300 * L (where 
L= 7.5 pe). This distance of travel by the shock front (~0.75 pz) is less 
than the diameter of the cluster cloud (2.0 Pc}yas we would expect. 


The effect of this impact layer is to convert high energy particles ints 
radiation in the visible and infra-red region at the spectrum. This, in turn, 
serves to increase the radiation density at these wavelengtns inside the cloud 
of the cluster. In orijer to determine if the cloud can handle this increase in 
radiation density at these wavelengths, we analyze the cooling capacity of just 


one mechanism - that of the dust grains. If we should find that an increase in 


radiation density would rsauire such 6 high temperature for the dust particles 


Cy 


that they could not remain in the sclid state, then our cloud could possibly 

be raised to a high equilibrium temperature. We pursue ‘this analysis in a way 
which seeks merely to set a limit as to the amount of cooling the dust can handle 
with respect to the increased energy input. 


From Roberts's data we have that the cloud absorbs about 5 magnitudes in 


the visible region of the spectrum (X~5x1L077 em). If we assume that the dust 


hes ad Lemperature of about 50 °K (ato er the cloud will be self-gravitating), 
then from Wien's law: » (em) = one) 
Max o rT Pie ) 


nN ( kK) ~ 
max (at 50 ee ec Ou acem 
since the absorption law for small (~107 cm) dust particles shows that absorp- 


tion varies inversely as wavelength, we conclude that effective absorption 


ee 








0% 


5 
ey DO “K is approximately 5 mag| = “20.05 magnitude. This means that 





LO 
the cloud has an absorption coefficient, (k)ae 5%, at the relevant wave- 


length, » | | . 


° 


Therefore, we apply the Stefan-Boltzmann radiation law (with the additional 
factor k) to our cloud, and thereby determine how much energy per second the 
Gleb in the cloud can radiate at 50 OK 


L(ergs/sec) = har-o Toi 


: -2,. 4 al 
where r=1pc, @= 5. (KU “ere cm deg sec ', P 
= Soe 
‘T = ao K anc kK = a) aie 
Piestituting these valves, we fiad that Le 2x10 % ergs ses (cooling rate that 


e220 handled ov a Roberts oloud with dust part 
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Now if we ccmpare this with the maxinnm possible energy ee from the cobpicien 
process, 6.4x10°~ ergs/sec, we conclude that the cloud can absorb this addi- 
fmemal Cnergy 1npul without vaporizing the Cus To 

By comparing with the above, the energy input expecced 1m the collision er 


~ 


va 
eeeiuster cloud with the intercloud medium of O.1 om 7 concentration veecene 


= 


fetes tnat the Robeyvts cloud will not be heated appreciably SxsSspLt for a cain 
Perea the collision junzturs. 

We also find that the IdLis cloud can withstand the thermal effects of 
cloud impact in the galactic disk. 

By a Similar analysis as used for the 260 © cloud, we find that the 


2 in the Idlis cloud, making D2Ga0- ae 


proton concentration is also loner 
Using the expression for transfer of momentum per unit area, we find that the 
a 
1onized layer travels at 0.3 km/sec through the globule (oy ot (500)"-10 
globule 10° 
Since this is below the speed of sound (a = 1.2 km/sec for the material), 


we ao not have a shock phenomenon in this case. 


ae 











CHAPTER V 


DISCUSSION OF RESULTS AND CONCLUSIONS 


From the calculations in Chapter IV, we can conclude that a globular cluster 
will not necessarily lose its interstellar material in passing through the galactic 
disk. 

We have seen that, for a cloud of the proportions described by Roberts and 
aiso for the smaller, denser globule of Idlis, momentum transfer will be sufficient 
in about half the cases to eject it from the cluster. We also learned that a 
shock front forms during collision with clouds of the galactic disk and that the 
recombination of hydrogen, and. redia ~icnm in the inmire-red by dust gcreaiis 18 a0ce.- 
ficient cooling process to dissipate the energy of collision which, therefore, 
will not cause heating and expansion of the cloud. 

From considerations discussed in Chapter III, we cunclided that the center 


of the cluster was the most satisfactor 


a 
T 
2 
eae 
bs 
O 
t~ 


ronment for HI regions. 


a © 


Accordingly, it would seem plausible that dirfuse material accumulates in 
Mam Center of the cluster and is excited into orbits at various distances on 
passages through the galactic disk. We would not expect to see those clouds in 
toe cluster at great distances from the center, even if they occurred there, nor 
those near the center on photographs with exposures typical of those used by 
Roberts. 

We have examined at some length possible correlations for those 12 clusters 
(which show obscuring regions) out of the 32 examined by Roberts. Correlations 
could be expected for at least two reasons: (1) differences in interstellar 
-avironments and evolutionary activity among clusters, and (2) significant, 
-ystemstic differences in the structure of the galactic disk. We must use caution 
4 attempting a correlation since it is not at all certain that the 20 clusters 
Without reported obscuration contain less diffuse material. Roberts has stated 


(private communication, 1963) that he was conservative in his estimates of 


oe 





[ee 











clouds and that his list of 12 could possibly be increased by 50%. 

We have noted that only about half of the globular clusters have been re- 
ported as having RR Lyrae variable stars (Hogg, 1959). Adding w Certauri 
(Fitzgerald 1955) and M4 (Idlis) to Roberts's list of 12 clisters, we note 
that RR Lyrae stars have been sighted in all but two. Of the 20 clusters 
ixamined by Roberts with negative results for obscuring clouds, 13 are listed 
os naving RR Lyrae stars. Generally speaking, the Lack of RR Lyrae stars in a 
cluster is attributed to the absence of stars in the appropriate evolutionary 
stage. While the above rigures don't indizate a significant dirference, in 
this case, some similar criterion having evolutionary Significance might possibly 

«+ uncovered by a closer search and classitication of variable stars and extremely 
blue stars in globulur clusters. 

We learned in Chapter LV that collisions with the dense clouds of the disk 
meeedeemajor factor for loss of clouds from the globular clusters. The force 
meroea Oy the intercloud méediue is net significart. Accordingly, we might 
expect that if there exists a large, contiguous area or tne galactic plane with 
& significantly smaller proportion ot dense (10 ca") structure than in the 
solar neighborhood, that those clusters passing througn this region in their 
galactic orbits would retain more interstellar material and be more likely to 
show obscurations. 

Rougoor and Oort (1960) have remarked that there appears to be a relatively 
small amount of dense HI structure in the region 600 pe to 3 kpc distance from 
the center in the galactic plane. Optical observations of the Sb galaxy 
icc 2841 (Sandage, 1961) show an amorphous region devoid of dust or spiral 
+ructure in the same region (as scaled to our Galaxy). 

We have determined, by the use of a formula from Eggen, Lynden-Bell, and 
Sandage (1962), that the average orbit of a globular cluster (e = 0.8) with 


° -O. ‘ : 
apogalacticon, 1O kpc, sweeps out about 205” between successive apogalacticon. 


259° 











we also have determined that our average 0.8 e orbit will pass within 3 parsecs 
or the center of the Galaxy during disk transit if the semi-major axis Lies 
within a cone with vertex angle of 50 eentered on the galactic axis. We 
have attempted a correlation, without success, tetween the present Location 
of clusters within the cone and the presence of obscuring material, The fact 
that heliocentric distances of globular clusters have probable errors up to 50% 
mekes any such correlation doubtful. 

Perhaps the two most significant recommendstions for the future study of 
globular clusters are: (1) An coservetionel program that will devermire vetter 
The Mamber of hot, blue stars in glotvilar clusters. For example, the incom- 


eteness factor of avout 15 ror 5 stars is. MS could ce improved by more care-= 


Sul observations. These stars are significazt in their ¢rrect on the lonization 


Ot HJ clouds, and the radiation field of these stars migkt be the most important 


Sweeresantiuencing the litetime of FI clowis, since they have a reasonamile chance 
“ Sirviving passage through the galactic disk. (2) It is desirable to have 


amore detailed program for observing obscured regions in glolilar clusters. 
Of primary interest would be the photographing of clusters at many dirferent 


exposures to determine better the cloud Locations and their absorption. 
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